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Saturation of gluon distribution is a consequence of the non-linear evo¬ 
lution equations of QCD. Saturation implies the existence of so called satu¬ 
ration momentum which is defined as a gluon density per unit rapidity per 
transverse area. At large energies for certain kinematical domains satura¬ 
tion momentum is the only scale for physical processes. As a consequence 
different observables exhibit geometrical scaling (GS). We discuss a number 
of examples of GS in different reactions. 
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1. Introduction 

At the eQCD meeting in 2013 jjj we have discussed the emergence of 
geometrical scaling [2j for F 2 /Q 2 in deep inelastic scattering (DIS) [3] and 
for charged particle distributions in proton collisions [3|. Here, after short 
reminder, we extend this analysis to (p r) (N c h) correlation [5j [6] and to 
heavy ion collisions (HI) [7j. References m , m-m include a more complete 
bibliography of the subject. 

Geometrical scaling hypothesis means that some observable a that in 
principle depends on two independent kinematical variables, say x and Q 2 , 
in fact depends only on a specific combination of them denoted as r: 

a(x,Q 2 ) = S ± F(r). (1) 

Here function F in Eq. 0 is a dimensionless function of scaling variable 

T = Q 2 /Ql{x). (2) 
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and 

Qs( x ) = Ql(x/x 0 )~ X (3) 

is the saturation scale. S± is a transverse area that corresponds to the 
overlap of hadrons colliding at fixed impact parameter b (or integrated over 
db), or - like in the case of DIS - it is a cross section for large dipole 
scattering on a proton. Qo and xq in Eq. © are free parameters, which can 
be extracted from the data within some specific model for a, and parameter 
A is a dynamical quantity of the order of A ~ 0.3. Here we shall test the 
hypothesis whether different pieces of data can be described by formula 
0 with constant A, and what is the range of transverse momenta where 
GS is working satisfactorily. Throughout this paper we shall be neglecting 
logarithmic energy dependence due to the running of a s . 


2. Deep Inelastic Scattering at HERA 

Let us start with DIS where the relevant scaling observable is F 2 (x)/Q 2 [2j. 
In Fig. [l] we plot F 2 (x)/Q 2 as a function of Q 2 (left panel) and in terms 
of r for A = 0.329 (right panel) for combined HERA data [8]. Different 
points correspond to different Bjorken x's. We see from Fig. [l] that points 
of different Bjorken x’s scale very well with some exception in the right part 
of Fig. 0b. These points, however, correspond to large Bjorken x’s where 
GS is supposed to break. 




Fig. 1. Combined DIS data [5] for F 2 /Q 2 . Different points forming a wide band as 
a function of Q 2 in the left panel correspond to different Bjorken x’s. They fall on 
a universal curve when plotted in terms of r (right panel). (Figure from the first 
paper of Ref. 0)- 






praszalowicz'GS'eQCD printed on May 12, 2015 


3 


3. Inelastic p t spectra at the LHC 

In hadronic collisions at c.m. energy W = \/s particles are produced in 
the scattering process of two patrons characterized by Bjorken x’s 


X\2 = e ±y pt/W. 


(4) 


For central rapidities x = x\ ~ X 2 - Geometrical scaling in this case means 
simply that [HJ: 


dN 


dyd 2 pT 


= S±F(t) 


(5) 


y~0 


where F is a universal dimensionless function of the scaling variable 


r = p\/Ql(x) = p\/Ql (pt/(x 0 W)) X . (6) 


In Fig. [2] we plot ALICE pp data PJ in terms of pt (left panel) and in 
terms of scaling variable r (right panel) for A = 0.22. We have found by 
a model independent analysis that the optimal exponent A = 0.22 — 0.24 
B33> which is smaller than in the case of DIS. Why this so, remains to be 
understood. 


,i=o , 1 = 0.22 




Fig. 2. Data for pp scattering from ALICE [9] plotted in terms of pt and -y/r. Full 
(black) circles correspond to IF = 7 TeV, down (red) triangles to 2.76 TeV and up 
(blue) triangles to 0.9 TeV. 
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An immediate consequence of GS for the px spectra is a power-like 
growth of multiplicity with energy. Indeed, since 

Pt = Q s (W)t ^ 2+a ) (7) 

where the average saturation scale is defined as 

Q B (W) = Qo (x 0 W/Q 0 ) x/{2+x) (8) 


one arrives at 


dN 

dy 


S ± Q s 2 (IG) x 


1 

2 +A 


Jt(t)t 2 ^ +x ^ 


(9) 


Data indeed support the power-like growth of inelastic multiplicity as s 01 
as predicted by GS by Eq. Q for A = 0.22 — 0.24. 


4. Mean pt in hadronic collisions at the LHC 



Fig. 3. Mean (j>t) in pp collisions at 7 TeV (full black circles), 2.76 TeV 
(full red down-triangles), 0.9 TeV (full blue up-triangles) and in pPb colli¬ 
sions at 5.02 TeV (full brown diamonds) plotted in terms of scaling variable 
(W/W 0 ) x ^ 2+x ^N c l/S ± . For pp Wo = 7 TeV and for pPb W 0 = 5.02 TeV. 
(Figure from the second paper of Ref. [5].) 

Another consequence of Eq. © is that 0 

(pt) ~ Qs(W) , (10) 

which means that (pp) rises with energy as W x ^ 2+x \ which is indeed seen 
in the data. On the other hand, since the saturation momentum is by 
Eq. ([8]) equal to the gluon density per transverse area, one easily derive 
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the correlation between mean py and charged particles multiplicity at given 
energy W 


W \ A /( 2+A ) 

(FT)\ W ~ ( Wa ) 


N A 


*-bL (-^Vch ) | \y 0 


( 11 ) 


By fixing multiplicity, one is probing some fixed impact parameter that 
corresponds to the overlap transverse area S± (N c h ) that itself is by con¬ 
struction both multiplicity and energy dependent. Therefore one needs a 
model for S±(N C h). To this end we have used the Color Glass Condensate 
result for pp and pA collisions m- The result is plotted in Fig. [3] where we 
plot ALICE data [T2j as a function of scaling variable defined in Eq. (11) 


5. Geometrical Scaling in heavy ion collisions 



Fig. 4. Illustration of geometrical scaling in heavy ion collisions at different energies 
and different centrality classes. Left panel shows charged particle distributions from 
ALICE [15], STAR jTGl |T7j and PHENIX (T8] [19] plotted as functions of fry ■ In 


the right panel the same distributions are scaled according to Eq. (131. 


GS for particle spectra in HI collisions has been already discussed in 
Ref. [7] and in Ref. m for photons. HI data are divided into centrality 
classes that select events within certain range of impact parameter b. In this 
case both transverse area S± and the saturation scale Ql acquire additional 
dependence on centrality that is characterized by an average number of 
participants IY part . We have |l3l fI3]: 


S± ~ and „ A p a r t • 


(12) 
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Therefore in HI collisions 


1 dN ch 
Nevt ^Vpart d^d 2 p T 


-jFfr) where r 

Q o 


Pt PPt\ x 
at 1 / 3 n 2 V W J 

1 v part Vo 


(13) 


In Fig. [4] we plot LHC and RHIC data in terms of px (left panel) and sjr 
for A = 0.3 (right panel). One can see an approximate scaling of, however, 
worse quality than in the pp case. 


To summarize: a wealth of data in hadronic collisions exhibit GS. This 
may be interpreted as a signature of saturation. However some details, like 
the non-universality of the value of A, remain to be understood. 

This work was supported by the Polish NCN grant 2014/13/B/ST2/02486. 
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